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ABSTRACT 

A three-dimensional radiative-dynamical-chemical model 
has been developed and used to study some aspects of mod- 
elling the polar lower stratosphere. The model includes 
a comprehensive gas-phase chemistry scheme as well as 
a treatment of heterogeneous reactions occurring on the 
surface of polar stratospheric clouds. Tracer transport is 
treated by an accurate, non-dispersive scheme with little 
diffusion suited to the representation of strong gradients. 
Results from a model simulation of early February 1990 
are presented and used to illustrate the importance of the 
model transport scheme, l he model simulation is also used 
to examine the potential for Arctic ozone destruction and 
the relative contributions of the chemical cycles responsi- 
ble. 

INTRODUCTION 

We present results from a 3D model which has been - de- 
veloped for the study of trace gases in the polar lower 
stratosphere. A detailed account oi the model is given in 
a paper in preparation (Chipperfield et al 1993) which de- 
scribes a series of experiments in the 1989-90 Arctic winter 
including a comparison with the available measurements. 
In this short article we concentrate on illustrating the im- 
proved model performance on changing from a spectral 
transport scheme to the second order moments scheme of 
Prather (1986) and also on the efficiency of various cat- 
alytic. ozone destruction cycles in the polar lower strato- 
sphere. 

The northern winter of 1989-90 was extremely cold with 
temperatures in the polar lower stratosphere below the sea- 
sonal average (e.g. see Naujukat et al 1990). An inspec- 
tion of the ECMWF (European Centre for Medium-Range 
Weather Forecasts) analyses for this winter shows that type 

I (nitric acid trihydrate) PSC’s could have formed from 
mid December to early February at which time a stiato- 
spheric wanning occurred ( although the voitex did not. 
break down). Temperatures in late .Januaiv and early 
February were cold enough to permit the foimahou of tv pc 

II (ice) PSC’s. The persistence of cold temperatures into 
February would have provided the greatest opportunity of 


PSC-processed air to experience sunlight, especially as dur- 
ing this period the cold vortex was displaced from the pole 
towards northern Scandinavia. 


THE MODEL 


We have used a middle atmosphere general circulation 
model (GCM ) to generate wind fields which are used in a 
chemical transport model (CTM). The GCM used is the 
tropospheric-stratospheric version of the 'Emeraude spec- 
tral model described by Cariolle et al (1990). The model 
has 30 levels (with a resolution of about 1 . 8 km in the lower 
stratosphere) and, for the experiments presented here, has 
been employed in the horizontal truncation of T21 (corre- 
sponding to a gaussian grid of 5.6° x 5.6°). The CTM uses 
the same grid as the GCM and advects the tracers using the 
scheme of Prather (1986) which conserves the second order 
moments of the tracer distribution. The chemical species 
are transported using this scheme rather than the spectral 
scheme of the GCM as it is less diffusive and better at 
representing the strong gradients often found in chemical 
species concentrations, for instance at the vortex bound- 
ary (see below). The CTM contains 20 integrated chem- 
ical species (O 3 , 0. H 2 O 2 , NO, NO 2 , NO 3 , N 2 O 5 , HNO 3 , 
HO 2 NO 2 , Cl, CTO. CION O- 2 , HOCT. HC1, OCR), Cl a 0 2 , 
Br, BrO. BrONO>. Hi C l) and 3 which are assumed to be 
in photochemical equilibrium (H. OH. HO 2 J. Photochem- 
ical data is taken from NASA/JPL (1990). The chemical 
integration is performed using the semi-implicit symmet- 
ric scheme of Ramaroson et al (1992). hixed zonal mean 
fields of CH 4 and CO (taken from a 2D model) are used 
for reactions such as Cl + CH.j — + HO + CH 3 . The model 
temperature, water vapour and nitric acid fields are used 
to predict the occurrence of type 1 (Hanson and Mauers- 
berger 1988) and type II PSC’s (Murray 1967). When such 
PSC’s occur the following five heterogeneous reactions are 
activated: 


N 205 (g) + H^Ofs) 
NUO- (g) t HCl(s) 
( TON 0 2 (g) -r H 2 0 ( s ) 
CIONO 2(g) t HCliVi 


-► 2HN0 3 (s) 

- (TNOjfg) -r HNO :3 (.s) 

— HOCl(g) + HNOj(s) 
— > CCig) + HNO;j(s) 


496 

C'X 


y 



HOCl(g) + HCl(s) -> Cl a (g) + HN0 3 (s) 

The rates of these heterogeneous reactions are calculated 
using reaction probabilities from experimental data (see 
WMO/UNEP 1991 for discussion and recommended val- 
ues). 

INITIALISATION 

The GCM is initialised with meteorological data from the 
ECMWF below lOhPa and from climatology above. The 
0 3 field was generated from TOMS data using the method 
described by Riishojgaard et al (1992). The other chemical 
species were initialised using data from a 2D model. Pole- 
wards of 30°N, between 200hPa and lOhPa, the longer-lived 
chemical species were initialised with the 2D data trans- 
formed into the coordinates of potential vorticity (PV) and 
potential temperature (0) in order to capture the azonal 
nature of the polar vortex. 

RESULTS 

To illustrate some aspects of the model performance we 
present results from four 10 day simulations of the CTM 
which start on the 4th February 1990. During this period 
temperatures in the polar vortex were cold enough to al- 
low the formation of both type I and type II PSC ! s and 
the vortex itself was displaced towards Scandinavia (figure 
1). The winds for the CTM were obtained from a 10 day 
run of the GCM starting from the same date. The low res- 
olution (T21) of the GCM run means that after a few days 
the meteorological forecast of the model is degraded. Run 
A of the CTM contained only gas-phase chemistry while 
run B included the treatment of heterogeneous reactions 
described above. To compare the effect of two transport 
schemes we performed two further runs (C and D) of the 
CTM (equivalent to runs A and B respectively) using the 
spectral scheme (as employed in the GCM) for tracer ad- 


vection. In runs C and D negative mixing ratios caused by 
the transport scheme were corrected by redistributing the 
negative mass amongst other members of the same chemi- 
cal family within a model box. 

Chlorine Processing 

Figure 2a shows the difference in HC1 mixing ratio on 
the 500K isentropic surface between run B and run A after 
10 days of model integration. Air within the polar vortex 
(indicated by the PV map in figure 3) has been processed 
by the localised region of PSC’s over northern Scandinavia 
in run B, converting virtually all of the HC1 (2ppbv) at this 
altitude to active chlorine (the initial conditions assumed 
that the chlorine species were partitioned purely on the 
basis of gas-phase chemistry). Note the sharp gradients in 
HC1 at the edge of the PSC-processed region in figure 2a. 
The non-diffusive nature of the Prather transport scheme 
performs well in maintaining these strong gradients. Fig- 
ure 2b shows the difference in mixing ratio for C10N0 2 
At the centre of the vortex is a minimum surrounded by 
a ring of high C10N0 2 . In calculating the rates of chemi- 
cal reactions the model simply uses the average concentra- 
tion of each tracer within a box, which therefore assumes 
rapid, complete mixing within that box. Therefore in run 
B the rate of production of C10N0 2 is probably overes- 
timated. To illustrate the improvement of this transport 
scheme over the spectral scheme, as employed in the GCM, 
figure 4a is the same as figure 2a but for runs C and D. 
Although the processing of HC1 in the vortex is similar the 
gradients at the edge of the vortex are smoothed out and, 
more unrealistically, the HCl-poor air is Transported’ to 
mid- latitudes. The limitations of the spectral scheme, in 
the comparatively low resolution of T21, are more apparent 
with ClONOj. Figure 4b shows much larger mixing ratios 
of C10N0 2 than figure 2b both inside and outside the po- 
lar vortex. A local minimum in the difference is coincident 
with the region of PSC’s at 8CPN. 90° E. Chlorine nitrate 



Fig. i. Temperature (K) at oOhPa on the fit !i f ebruary 
1990 from the T21 GCM simulation (initialised on Febru- 
ary 4th) used tu produce tin* wind fields for the CTM. At 
this altitude type 1 PSC s could torm at 19 )K and type II 
PSC’s at 190K. Contour interval 2K. 


Fig. 2. Difference in a) HC1 and b) C10X0 2 mixing ratio 
( ppbv) between run B and run A on the 50QK surface on 
14lh February. Contour interval 0.2pphv. 
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Fig. 3. Potential vorticity (xlO 6 Kiras'" 1 kg" 1 ) on the 500K 
surface on 14th February 1990 from the GCM model sim- 
ulation. Contour interval 10xl0 6 . 


Fig. 4. Difference in a) HC1 and b) C10N0 2 mixing ratio 
(ppbv) between run D and run C on the 500K surface on 
14th February. Contour interval 0.2ppbv. 


is produced by reaction of CIO (typically high inside the 
vortex) with N0 2 (typically high at mid-latitudes). Figure 
4b shows the dangers of the spectral transport schemes 
where diffusion, due to the Gibb’s effect, ‘transports’ CIO 
to mid-latitudes and N0 2 into the polar vortex resulting 
in the formation of large amounts of C10N0 2 . The prop- 
erties of various transport schemes have been well investi- 
gated in several ID passive tracer experiments (e.g. Muller 
1992 and references therein). However, the examples above 
show that when chemical interactions between tracers are 
considered the problems are compounded and results from 
a diffusive scheme, such as the spectral one employed here, 
rapidly become unrealistic. The performance of the spec- 
tral scheme could be improved by the use of a more sophis- 
ticated filter to remove the Gibb’s fringes, or by increasing 
the resolution. 


Ozone budget 

The processing of air by PSC’s in run B leads to elevated 
levels of CIO and BrO in the polar vortex. Figure 5a is a 
latitude height cross section of CIO at 33°E at 12H GMT 
on February 6th. A peak mixing ratio of l.Oppbv of CIO 
is located at 30hPa between 55° N and 70°N. The altitude 
distribution of the CIO is determined by the upper limit of 
PSC’s in the model (around 20hPa) and the abundance of 
inorganic chlorine at lower altitudes. In the Arctic, clima- 




Fig. 5. Mixing ratio of a) CIO (ppbv) and b) BrO (pptv) 
at 33° E on 6th February at 12H CM I from run B. 


tology shows that PSC’s tend to occur at higher altitudes 
than in the Antarctic. The peak mixing ratio of BrO (figure 
5b) is Spptv at, a similar altitude and between 50°N and 
65°N. The subsequent destruction of ozone in. the model 
has been analysed for the following 3 cycles, which domi- 
nate the ozone loss in the polar lower stratosphere: 


I 


CIO + CIO + M 
C1 2 0 2 T hu 
C10 2 + M 
2( Cl T 0 3 
net : 2O3 


C1 2 0 2 + M 
-+ C10 2 + C1 

- ci + o 2 + m 

- C10 + 0 2 ) 

- 30 2 


II 


CIO + BrO -4 
CI + O3 -4 
Br + 0 3 — ► 
net : 2O3 — ► 


C1 + Br + 0 2 
CIO 4* 0 2 
BrO + 0 2 
30 2 


III CIO + O 

CI + O3 
net : 0 T O3 


Cl + 0 2 
C10 + 0 2 
20 2 


Cycle II could also be initialised by the formation and 
subsequent photolysis of BrCl. The instantaneous rate of 
ozone loss due to these 3 cycles is shown in figure 6. The 
efficiency of these cycles depends strongly on altitude and 
also on the latitude. Destruction due to cycle I (figure 6a) 
peaks at over 25ppbv/day at 40hPa, 60° X. The efficiency of 
this cycle is greater at lower altitudes due to the 3 body for- 
mation of C1 2 0 2 . The low abundance of CIO below 50hPa, 
however, limits ozone destruction due to this cycle below 
this height. Loss due to cycle II is similar to that due 
to cycle I and shows a maximum of over 20ppbv/day at 
around 40hPa. Ozone loss at lower altitudes would have 
a larger effect on the integrated column but on the model 
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Fig. 6. Ozone loss (ppbv/day) due to a) cycle I b) cycle II c) cycle III and d) 
production due to photolysis of O 2 at 33°E on 6th February at 12H GMT from 
run B (note irregular contour spacing in figures 6c and 6d). 


day shown the catalytic cycle with the largest rate of ozone 
mixing ratio loss is cycle III. Between 20-30hPa loss due to 
this cycle is over 35ppbv/day. The efficiency of cycle III in- 
creases with height as the concentration of atomic oxygen 
increases and in figure 6 it is the dominant loss mechanism 
above 30hPa (in the trajectory model study of McConnell 
et al (1991) cycle III was the major loss cycle above 19km). 
To some extent loss due to cycle III at around 20hPa is off- 
set by in-situ production of odd oxygen by the photolysis 
of O 2 (figure 6d) which increases rapidly at higher altitudes 
and lower latitudes. At 20hPa, 60°N the rate of production 
is 50ppbv/day. An estimate of the total loss of ozone dur- 
ing the winter period would need to take to account of the 
efficiency of these catalytic cycles as a function of time and 
location. This is discussed in more detail in Chipperfield 
et al (1993). 

SUMMARY 

These results highlight the importance of a non-diffusive 
transport scheme in 3D chemical model studies of the po- 
lar stratosphere. The expense of 3D chemical calculations 
means that comparatively low resolutions (such as T21) are 
usually employed. Such low resolutions provide a stern test 
for the transport scheme employed especially when chemi- 
cal interactions between species are considered. In our ex- 
periments the Prather (1986) transport scheme performs 
significantly better than a spectral scheme. The perfor- 
mance of the Prather scheme itself could be improved by 
taking account of the species distribution within a box in 
the calculation of chemical rates. The chemical cycles con- 
trolling the destruction of ozone in the Arctic lower strato- 
sphere vary strongly with altitude. The rate of ozone loss 
due to cycles involving CIO + BrO can be significant com- 
pared to loss due to CIO + CIO. When PSC processing 
leads to elevated CIO at altitudes above about 30hPa the 
destruction cycle involving CIO + 0 becomes efficient. 
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